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ABSTRACT 

In previous work, we discussed the detection of reverberation delays in the NLS1 
1H0707-495. The delays originate close to the black hole. Here, we show that 
RE J1034+396 shows very similar lag properties. At low frequencies (< 1 x 10 -4 
Hz), the time lag between energy bands increases with energy separation, similar to 
that commonly seen in Galactic black holes and other AGN. At higher frequencies 
(~ 3.5 x 10" 4 Hz), the soft (< 1 keV) and hard (> 3 keV) bands lag behind the 
intermediate band (1-3 keV). The simplest interpretation is that the intermediate 
band is dominated by the direct power-law continuum, while the soft and hard bands 
are dominated by the relativistically-smeared reflected emission. The low frequency 
delays are present in both available observations. The high frequency lags are only 
seen in one observation. In the observation where high frequency reverberation delays 
are observed, the spectrum contains a power-law component and there is a QPO in 
the light curve. In the other observation, no power-law component is required and no 
QPO is seen. The lags originate a few gravitational radii from the black hole, and the 
QPO is associated with the power-law emitting corona. 

Key words: X-rays: galaxies - galaxies: individual: RE J1034+396 - galaxies: active 
- galaxies: Seyfert - galaxies: nuclei. 



1 INTRODUCTION 

RE J1034+396 is a Narrow Line Seyfert 1 galaxy (NLS1, 
z = 0.042). It has an unusual multi- wavelength spectrum, 
with no significant big blue bump in the UV/Optical bands. 
It is very bright in the EUV and soft X-rays making it 
one of the brightest objects with a so-called 'soft excess' 
(|Puchnarewicz et al.|[l995l : IPounds et alJll995h . The soft ex- 
cess below rsj 2 keV is a common feature in the X-ray spectra 
of many AGN, particularly of the NLS1 class. It is too hot to 
be thermal blackbody emission from the accretion disc, while 
the similarity of its shape in objects with different masses 
argues against an origin in a cool Comptonised region exist- 
ing along with the hot corona emitting at higher energies. 
A corona of this type is expected to depend on the seed 
photons from t he disc, which in turn dep end on the mass of 
the black hole ([Gierlinski fe Pond 12004 ). Atomic processes 
are more likely to be at work through line-of-sight absorp- 
tion or reflection. The smoothness of the ob served emission 
rules out the former ([Schurch &; Donell2008l ). while the lat- 
ter works if the emission originates very close to the black 
hole where stron g relativistic effects acts to smear ou t sharp 
atomic features ([Crummy et al.l 2006: iNardini et aj]|201lh . 
One other object, 1H0707-495, also shows a strong 
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soft excess below 1 keV. It was shown recently that 
because of the high iron abundance, its s oft excess is 
main l y due to a stron g iron L emission ([Fabian et al.l 
2009; Zoghbi et al. 201 0|) produced by partially ionis ed re- 
flect ion^^^^^bia^lS^; [Georg^ The 
compactness of the reflecting region and its proximity to 
the black hole produces strong relativistic broadening that 
smears out the reflection features. The iron L-line is ob- 
served t ogether with the more commonly seen broad K line 
(e.g. see Miller 2007 for a review). 

1H0707-495 also shows a time delay between the hard 
band (1.0-4.0 keV) dominated by the direct power-law emis- 
sion, and the soft band (0.3-1.0 keV) dominated by th e 

;IZq£ 



reflected emission ([Fabian et al.l 120091 : IZoghbi et al.1 l2010h . 
The soft lag is the first to be detected, and is a confirmation 
that the soft band is dominated by reprocessed emission. 
The magnitude of the lag itself (30-50 seconds) is consis- 
tent with light travel time at ~ 2 gravitational radii from 
the black hole, as inferred from the energy spectrum. De- 
tailed energy-dependent lag spectra show that the lag traces 
the shape of t he reflection spectru m, further confirming its 
interpretation ([Zoghbi et al.ll201ll ). 

In this work, we investigate the nature of the soft excess 
in RE J1034+396 drawing on the similarities with 1H0707- 
495. We mainly concentrate on time lags and their inter- 
pretation. Spectral fitting and spectral variability properties 
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have already been investigated in iMiddleton et al.1 (|2009h . In 
that work, the increase of rms variability with energy was 
the basis of favouring a low-temperature Comptonisation 
interpretation for the soft excess. This component remains 
constant while the variability is produced by the hard tail 
changing in normalisation. Although such low-temperature 
Comptonisation is not generally seen in AGN, they suggest 
that RE J1034+396 is accreting at super-Eddington rates. 
It is therefore exceptional, which could also explain the fact 
that it is the only AGN to show quasi-periodic oscillations 
in its light curve. 



2 OBSERVATIONS & DATA REDUCTION 

We use two publicly available XMM- Newton observations 
of RE J1034+396. These were taken in 2007 (obs. id 
0506440101) and 2009 (obs. id 0561580201). The EPIC-PN 
detector was used in full and small window modes in the 
first and seconds observations respectively. The total expo- 
sure after removing high background periods was about 80 
and 50 ks for the first and second observations respectively. 
The first observation was piled- up. The spectra were ex- 
tracted in different ways to test the effect of pileup. First, 
we used PATTERN=0 events only, and although this may 
give the wrong total number of counts, the spectral shape 
is not distorted. Also, we extracted the spectra using dif- 
ferent extraction regions where the central highly piled-up 
parts of the Point Spread Function (PSF) were excluded. In 
the spectral fitting, spectra extracted in the different ways 
have minor differences. Although for the purpose of inter- 
preting the lags, we are only interested in a general spectral 
model, very good fits are found regardless of the type of ex- 
traction, but with parameters that differ slightly which is 
irrelevant for the this work. We also calculated the lag with 
and without removing the central regions of the PSF. These 
are discussed further in Sec. 13.11 



3 TIME LAGS 

3.1 Lag measurements 

It has been known that lags in black holes de pend on both 
frequ ency and energy at the same t ime (e.g. iNowak et al.l 
1999 for galactic black holes and McHar dvet al.1 20071 : 
Arevalo et all 120081 for AGN). In Zogh bi et al.1 (|201lh . we 



showed that the high frequency lags in 1H0707-495 (rever- 
beration lags) are closely linked to the energy spectrum. In 
order to measure any time delays, the energy bands have 
to be selected appropriately. For RE J 1034+396 however, 
the components of the energy spectrum are not well known. 
The whole energy band (0.3-10 keV) and frequency bands 
need to be explored. The two observations are so different 
in terms of their spectra and light curve properties, that we 
explore them separately. 

Light curves were extracted in 50 energy bands with 
time bins of 200 seconds. These were selected so that they 
have roughly equal number of counts. Fourier-resolved time 
lags were calculated between each light curve and a light 
curve constructed from the whole 0.3-10 keV band, exclud- 
ing the current band. Taking the whole band as a refer- 
ence instead of a single individual band maximises the sig- 
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Figure 1. 2D plot of the lag as a function of energy and frequency. 
The lags have been smoothed to help identify the general trends 
in the lag spectra. The plot is for the first observation only as 
the two observations have different lag properties (see text for 
details). 



nal to noise ratio, while removing the current band from the 
reference ensures t he noise remains uncorrelated betwee n 
the bands (e.g. see IZoghbi et al.1 120 '111 : lUttlev et al.1 120 111 ). 
Time lags and their errors were c alculated following stan- 
dard methods ([Nowak et al.lll999l ). The frequency axis was 
binned in log-steps. 

Fig. [T] shows a 2d plot of the lag as a function of en- 
ergy and frequency. The image has been smoothed, as its 
purpose is only to show the general trends in the energy- 
frequency parameter space. At low frequency (left), the lag 
tends to increase with energy. At intermediate frequencies 
(~ 3.5 x 10 -4 Hz), the lag is negative around 2 keV. Now 
given this rough idea of the energy- frequency dependence, 
we investigate it further. 

Fig. [2] shows the energy dependence of the lag for two 
frequencies (< 10" 4 and (0.2 - 1) X 10" 3 Hz). Light curves 
are extracted in six energy bands, defined to roughly have 
the same number of counts. The values of the lag have 
been shifted so that the minimum is at zero. Because the 
lag is calculated taking the entire energy band as a refer- 
ence, only the difference between the bands and the shape 
are of importance. The shapes are clearly different in the 
two frequencies, an d roughly resemble those of 1H0707-495 
(Zog hbi et al.ll201lh . At low frequencies, the lag increases 
with energy. This is t he behaviour seen in several Galactic 
accreting black holes ([Miyamoto &; Kitamotol [l989h . Below 
1 keV, the lag is more constant, indicating the presence of 
a separate spectral component. At higher frequencies, the 
lag dependence is reversed, and now decreases with energy, 
with a turn back at ~ 4 keV. 

The ener gy spectrum has at le ast two separate com- 
ponents (e.g. Mid dleton et al.ll20Q9h . dominating the spec- 
trum above and below 1 keV respectively. This appears to 
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Figure 2. Lag dependence on energy for the first observation of 
RE J1034+396 at high ((0.2-1) x lO" 3 Hz, top) and low (< 10~ 4 
Hz, bottom) frequencies. Lags are calculated with respected to the 
entire energy band (0.3-10 keV), then shifted so the minimum lag 
is zero. This shift implies the lag values in the top and bottom 
panels represents lags with respect to different bands. 



be consistent with the low frequency lags. More interestingly 
however, there is a soft lag between the soft (< 1 keV) and 
harder band (> 1 keV), in a sense that the hard band leads 
the variability at ~ 3.5 x 10 
to 1H0707-495 (|Zoghbi et al.ll201ll 
the soft band (< 1 keV) is dominated by reprocessed emis- 



4 Hz . This is strikingly similar 
and would indicate that 



sion, while the hard band is dominated by a direct power-law 
component. The spectrum is investigated in detail in Sec.[4l 



3.2 Lag significance 

In order to check the lag significance, we first explore the 
coherence between the bands at the frequency of interest. 
The coherence function is a measure of the linear correla- 
tion between any two b ands at the specified frequency (see 
IVaughan &; Nowakll 19971 ) . The soft lag is maximum between 
^0.5 and 1.5 — 2 keV, and we concentrate on measuring the 
coherence between these bands (0.4-0.6 and 1.5-2 keV). Fig. 
[3] shows the estimated lag and coherence between the 0.4- 
0.6 and 1.5-2.0 keV energy bands. The lag is positive at 
low frequencies (hard lags) and negative at high frequencies 
(soft lags). The corresponding coherence shows a significant 
coherence between the bands below ~ 8 x 10 -4 Hz. The co- 
herence appears to drop slightly at 10 -4 Hz. This is again 
similar to 1H0707-495, and could be related to a transition 
between two variability processes dominating the low and 
high frequencies. 

The frequency of the soft lags in RE J1034+396 (~ 
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Figure 3. Top: Lag as a function of frequency between 0.4- 
0.6 and 1.5-2.0 keV bands. Bottom: The corresponding esti- 
mated coherence function . The coherence was measured following 
IVaughan fe Nowaldll997l . 
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Figure 4. Lag as a function of energy for the first XMM obser- 
vation at (0.2 — 1) x 10 -3 Hz showing the effect of pile-up on 
the lag. Panels a,b,c correspond to the lags measured when the 
source region is an annulus with inner radii of 0, 7 and 32 arcsec- 
onds respectively. Panel (d) shows the same plot from a simulated 
piled-up event file using SCISIM where no lag is expected. 
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3.5 x 1CT 4 Hz) compar ed to those in 1H0707-495 (~ 2 x 1CT 3 
Hz; IZoghbi et al.ll201lh reduces the effect of Poisson noise on 
the lag signal. We checked this using Monte Carlo simula- 
tions of light curves that have the same rms, inter-band lag 
and coherence as those o bserved. We followed the method of 
iTimmer fc Koenid (Il99gv_ The details of the simulation are 
very similar to those in IZoghbi et afl fcoiO). The lags in the 
simulation are recovered and their RMS spread is consistent 
with the measured errors. 

The first observation is highly piled- up. We repeated 
the lag calculations for several extraction regions, excluding 
the central parts of the PSF (using annuli with inner radii of 
0,7,14 and 32 arcsec). This has very little effect on the shape 
of the lag-energy plot. The results are shown in Fig. [4j The 
plots are for annuli of inner radius of (a), 7 (b) and 32 (c) 
arcsec. To make a further test, we used the XMM- Newton 
SCISIM packages to simulate a whole piled-up observation. 
The input of the simulation is a point source with a flux 
of 8 x 1CT 12 ergs cm -2 s -1 in the 0.3-10 keV band and a 
spectrum taken as a two power-law component model fitted 
to the data. The simulated events file were reduced in the 
normal procedure (using SAS vlO.0.0). The pile-up was of 
the same strength as in the real observation. The lag-energy 
dependence of the simulated data is plotted in Fig. SJ-d. The 
lag is constant as expected, and pile- up does not effect the 
lag measured at the frequencies of interest. 
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Figure 5. The spectra from the two XMM observations. Top: 
Unfolded spectra produced in a model-independent way following 
Nowak et al (2005). The bottom panel shows the ratio of the 
spectra to a model consisting of a sum of two power-law models. 



4 ENERGY SPECTRUM 
4.1 Spectral fitting 

RE J 1034+396 has one of th e strongest soft excesses known. 
It was studied extensively in Middle ton et al.l (|2009h . Using 
constraints from the RMS spectra, a Comptonisation model 
was favoured. Here, we explore the spectra from both obser- 
vations, and include information from the lag spectra. 

Fig. [5] shows the spectra from the two observations. 
They are plotted as unfolde d spectra (top ) in a model inde- 
pendent way following iNowak et al.l (|2005h . This corrects for 
the detector effective area and shows the main features in 
the spectrum. The bottom plot shows the ratio of the data 
to a double power-law model. All parameters of the model 
have been allowed to vary in the fit. The aim of these plots 
is to highlight the spectral changes and differences between 
the two observations. The spectrum not only appears to be 
pivoting around the 1 keV point indicating separate compo- 
nents above and below 1 keV, but also seem to have struc- 
tural changes below 1 keV that could be related to atomic 
features in the spectrum of the soft band. This is possibly 
an indication that the soft component (below 1 keV) is not 
a smooth power- law- like component, but rather it has some 
atomic features and the difference between the observations 
is partly due to ionisation changes. 

Atomic features in the spectrum (or changes thereof) 
can be due to absorption, reflection or a combination of 
both. The smoothness of the spectrum itself requires the 
features to be smeared out by high velocity shear in a wind 
or by relativistic broadening. Smoothing ab sorption features 
requi res unrealistically high velocities (~ c JSchurch &; Done] 
2008), so a remaining possibility is relativistic broadening 
of reflection emission, which is also motivated by the lag 



resul ts. We used a model similar to that u sed for 1H0707- 
495 dFabian et al J 120091 : IZoghbi et aDl2010h consisting of a 
power- law and blurred reflection. We used th e ionised reflec- 
tion table model from iRoss fe Fabianll2005l convolved with 
relativistic kernel kdblur (|Laorl 199 lh . Galacti c column in 
the d irection of the source is 1.52 x 10 20 cm" 2 (iStark et al.l 
I1992L modelled in Xspec with tbabs. lWilms et alJl2000h . To 
obtain more constraints at low energies, we also included an 
earlier observation from the ROSAT PSPC archives that ex- 
tends down to 0.1 keV. This provides better constraints on 
the accretion disc thermal emission peaki ng in the EUV and 
expe cted to contribute to soft X-rays (|Puchnarewicz et"al] 
Il995h . 

For initial modelling, we fitted the XMM-Newton data 
with a model consisting of a power-law plus blurred reflec- 
tion. This simple model did not give a very good fit. The 
strongest residuals were an edge-like feature at ~ 0.5 keV. 
Adding an edge to the model improved the fit significantly 
with the edge having an energy of 0.49 keV. This corre- 
sponds to oxygen photoelectric absorption at the source 
redshift (z = 0.042). The fact that it is variable indicate a 
change in the line-of-sight column density at the host galaxy. 
To model this, we included a photoelectric absorption com- 
ponent with variable abundance (tbvarabs). The redshift 
was fixed at that of the host galaxy. This improved the fit 
significantly (A% 2 > 100 per extra degree of freedom). To 
obtain further improvements, we also added data from the 
ROSAT PSPC, so the model now extends down to 0.1 keV, 
which helps constrain the thermal emission from the accre- 
tion disc, and any contribution it may have in the XMM 
band. 

Fig. [Bl- (a) shows the spectrum from the two XMM obser- 
vation (black and red) and the ROSAT observation (green). 
The spectra are again unfolded in a model-independent man- 
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Figure 6. (a) The spectrum of RE J1034+396 from the two 
XMM observations (black and red), and from the ROSAT PSPC 
(green) unfolded following Nowak et al. (2005). The bottom panel 
is a plot of the fit residuals in terms of sigma. (b): The best fit 
model for the first XMM observation. The model consists of a 
low temperature blackbody, a power-law and a blurred reflection, 
(c): Similar to (b) for the second XMM observation. The colours 
are similar to those in Fig. [5] 

Parameter Value 

Tpower — law 2.65 zb 0.03 



1.90t°;£ (r 9 ) 



fij £2, £3 



6.24 ±0.5 

715l§|, 346l^,184tf (erg cm s" 1 ) 



Table 1. The fit parameters for the model shown in Fig. [6] ri n is 
the inner radius of the reflecting region, q is the emissivity index 
of the disc, modelled as oc r~ q . £ is the ionisation parameter. 



ner (jNowak et alJl20Q5h . The models are plotted on the right 
for the first (top) and second (bottom) XMM observations. 
The model provides a good fit to the data (% 2 = 698 for 611 
degrees of freedom) . The blackbody has little contribution in 
the XMM band and mainly contribute to energies less than 
0.3 keV, and it is in full agreemen t with multi- wavelength fit s 
using optical to soft X-rays data (|Puchnarewicz et al.l fl995). 
where there is no strong big blue bump in the optical/UV 
but the source has a strong emission in the EUV/Soft X-ray 
indicating hot accretion disc emission. The disc component 
in our case is mainly constrained by the ROSAT PSPC data, 
and remains constant between the three observations. The 
best fitting parameters for the whole model are presented in 
Table ED 



4.2 Modelling the lag 

The soft lag discussed in Sec. 13.11 should, when interepreted 
correctly, be related to the energy spectrum. The spectral 
model (Fig. [6J show that for the first observation, the rel- 
ative contribution of the reflection compared to the direct 
power-law initially decreases with energy, it reaches a mini- 
mum at rsj 2 keV, and increases again. This is precisely the 
behaviour of the lag in Fig. [2] To show this more clearly, we 
plot in Fig. [71 (top) the lag-energy dependence, along with a 
model representing the ratio of the reflection component to 
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Figure 7. Lag as a function of energy for the first (Top) and 
second observation (Bottom). The whole 0.3-10 keV band is used 
as a reference (removing the current band each time). The red 
line is a model constructed by taking the ratio of the power- 
law component in the best fit model to the total (power-law 4- 
reflection). The lines represent the best fit model. More energy 
bands are used here compared to Fig. [2] to enable better energy 
resolution. The energy bands are selected in equal steps in log 
space. 



the total spectrum (power- law + reflection) shown in Fig. [6] 
The constant offset was added to the model to account for 
the reference, and multiplied by a factor representing the in- 
trinsic lag between the two components. The latter is found 
to be 296 =b 80 seconds. 

The two XMM observations are different spectrally and 
in the light curve properties. The second observation is 
reflection-dominated, where the direct power-law is not ob- 
served. Several AGN have been observed sometimes to be 
in a reflection-dominate d state, particularly at low fluxes 
(e.g. Zo ghbi et al.ll2008l ). Although the second observation 
is shorter, we calculated the energy-dependent time lag and 
show it also in Fig. [3 (bottom). The data do not show the 
lags seen in the first observation and are consistent with 
a constant. This is exactly what is expected, as the direct 
emission is not seen, and the lag now is measured between 
different parts of the reflection component itself. Despite the 
high uncertainties in the second observation, the model fit- 
ted to the first can be ruled out in the second observation 
at around 95 percent confidence. 
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5 DISCUSSION 

We have shown that there is a soft lag between the hard 
and soft components in the light curves of RE J1034+396. 
The simplest interpretation is that the soft component (< 1 
keV) is reprocessed emission produced in a partially ionised 
reflector. This object is very similar to 1H0707-495 in its lag 
properties, and its spectrum can be well modelled with a 
reflection model dominating the soft band. The best fitting 
spectral model matches very well the lag-energy spectrum. 

The smoothness of the soft excess requires the sharp 
reflection features to be smeared out by strong relativistic 
effects, implying that the emission region is very close to the 
black hole. Spectral fitting indicate the inner radius of the 
emission region to be < 3 gravitational radii, with an emis- 
sivit y index of ^ 6, slig htly flatter than that of the 1H0707- 
495 feogfibi et aLl koiO). indicating a slightly extended emis- 
sion region. Thermal emission from the accretion disc also 
contribute to the soft X-ray band (0.1-0.4 keV). The X-ray 
and EUV flux is unusually higher than the flux in the optical 
and UV, indicating that the big blue bum p, due to the accre- 
tion disc, is shifted to higher energies (|Puchnarewicz et al.l 
Il995h , and a small contribution is expected in the XMM PN 
band. 

A low variability, low temperature Comptonised 
emission can reprodu ce the soft excess in this object 
(jMiddleton et alJl2009h . However, it would be difficult, in 
such a model, explaining the temperature of the seed pho- 
tons, which appears constant, independent of the black hole 
mass. It is also interesting to explore the long term spectral 
changes. The differences between the ROSAT and XMM 
data (|6]-a) indicate that the soft band is indeed dominated 
by broad spectral features that change with time, possibly 
due to geometry and/or ionisation changes. Furthermore, 
the detection of a soft lag where the soft band (0.5-1 keV) 
lags the medium band (1-4 keV), with an energy dependence 
that has a minimum around 2 keV, also strongly supports 
the interpretation of the soft emission being due to repro- 
cessing. 

It is interesting to explore the measured lag quanti- 
tatively too. The lag in 1H0707-495, after correcting for 
the contribution o f individual components, is ~ 70 seconds 
(jZoghbi et al.ll201lh at a temporal frequency of 1 — 2 x 10 3 
Hz. In RE J1034+396, the measured lag is ~ 290 seconds at 
a frequency of ^ 3 — 4 x 10 -4 Hz. There is a factor of ^ 4 
difference in the lag and the time-scale of the lag, indicat- 
ing possibly a factor of 4 difference in mass. The frequency 
at which the lag is measured corresponds to the viscous or 
thermal time-scales at the radius the emission, while the lag 
is a measure of the light travel time between the illuminat- 
ing source and the reflector. Although the exact numbers 
depend on the exact unknown geometry, the lags are consis- 
tent with mass estimates of 1 — 2 x 10 6 M and 4 — 8 x 10 6 
M for 1H0707-495 and RE J1034+396 respectively. 

RE J1034+396 is the only Seyfert galaxy known to 
show a quasi-periodic oscillation (QPO) in its light curve 
(jGierlihski et al.ll2008l : IMiddleton et alJl20Q9h . The QPO is 
present in the first XMM observation but not in the second. 
The spectra of the source in the two observations appear to 
be different (Fig. [6}. Our spectral fitting has shown that the 
spectrum in the second observation is reflection-dominated. 
In the first observation however, the fitting requires the pres- 



ence of a power-law component that contributes significantly 
around 2 keV. This indicates a possible relation between the 
QPO and the power-law component in the energy spectrum. 
This is also support ed by the shape of the rms spectrum of 
the QPO (Fig. 2 in IMiddleton et al J l2009h , where the rms 
powe r in the QPO peaks around 2 keV. iMaitra &; Miller! 
(2010) attributed the QPO to an orbiting cloud close to 
the black hole, producing mainly an absorption edge at ^ 8 
keV. However, the apparent edge change is more likely to 
be due to continuum change and not a change in the atomic 
features. The shape of the RMS spectrum rules out this 
possibility. Linking the QPO to the power-law component is 
in fact not surprising given what is known about QPOs in 
Galactic black holes (GBH), where the QPO is associated 
with the hard spectral tail and not the therm al disc compo- 
nent (e.g. see iRemillard &; McClintockll2006h . This however, 
does not necessarily mean that observing a power-law com- 
ponent implies detecting a QPO. 

QPOs in GBH come generally in two types: low- and 
high-frequency QPO s and are observed in pairs wi th a fre- 
quency ratio of 3:2 ( IRemillard &; McCrintock]|2006l and ref- 
erences therein). Identifying the equivalent of the QPO in 
RE J1034+396 is not clear, but taking its mass to be a few 
xlO 6 M0, makes the QPO more likel y to be the equivalent 
to the high-frequency type fHFQPO. fGierlihski etai1l2008l : 
IMiddleton et al.|[2QQ9h . A simple relation between mass and 
the fundamental frequency was suggested for GBH of the 
form: M — ty&lf^ M^ where fp is the fundamen tal fre- 
quency of the QPO (Remillard & McClintock 2006). If we 
assume it applies also to AGN, the mass of RE J1034+396 
would be in the range 6 — 10 x 10 6 Mq (depending on which 
harmonics the observations corresponds to), which is very 
consistent with the estimates from the lag. 

The lag in RE J 1034 extends over a wide frequency 
range, similar to 1H0707-495. The frequency of the QPO 
(a narrow frequency band by definition) is smaller than the 
frequency where the lag is maximum (Fig. [T]and[3]). The 
accretion disc is reflecting the random stochastic signal of 
the illuminating source, which is this case happens have a 
QPO, but generally it does not. The difference in frequency 
between the QPO and the lag is possibly an indication that 
the corona producing the QPO is slightly extended, while 
the compact reflection from the disc is emitted in a very 
small region and dominated by light-bending effects. This 
can also explain the absence of the QPO in the reflection 
component, only the inner parts of the variable power- law 
illuminate the inner regions of the disc and not the regions 
producing the QPO. 
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